NEWTONIAN DYNAMICAL SYSTEMS 
ADMITTING NORMAL BLOW-UP OF POINTS. 



RUSLAN A. SHARIPOV 

Abstract. Class of Newtonian dynamical systems admitting normal blow-up of 
points in Riemannian manifolds is considered. Geometric interpretation for weak 
normality condition, which arose earlier in the theory of dynamical systems admitting 
the normal shift of hypersurfaces, is found. 



1. Introduction. 

Let S be a hypersurface in Riemannian manifold M. One of the ways for de- 
forming S consists in shifting points, which constitute S, along trajectories of some 
Newtonian dynamical system. Such situation arises in describing the propagation 
of electromagnetic wave (light) in non-homogeneous media in the limit of geomet- 
ric optics. Hypersurface S models wave front set (the set of points with constant 
phase), while trajectories of shift model light beams. Newtonian dynamics of points 
of Riemannian manifold M in local coordinates x 1 , . . . , x n in M is described by a 
system of n ordinary differential equations 

n n 

(1.1) x k +J2'£r^x i ^ = F k (x\... ,x n ,x\... ,x n ), 

i=l j=l 

where k = 1, ... , n. Here T k j are components of metric connection for basic 
metric g of the manifold M. Quantities F k are components of force vector F. 
They determine force field of dynamical system (1.1). In the equations (1.1) they 
play role of perturbing factor, due to them trajectories of dynamical system (1.1) 
differ from that of geodesic flow for the metric g. 

At each point p of hypersurface S we fix some vector of initial velocity v(p) and 
determine trajectories coming out from all points of hypersurface S by setting the 
following Cauchy problem for the equations (1.1): 

(1.2) x k \ = x k (p), x k \ =v k (p). 

lt=0 lt=0 

Here v k {p) are the components of vector v(p) in local coordinates x , ... , x n . 
Having displaced for time t along trajectories determined by initial data (1.2), 
points of the hypersurface S constitute another hypersurface St- As a result we 
obtain a family of hypersurfaces and a family of shift maps ft: S — *• St, which are 
local diffeomorphisms for sufficiently small values of parameter t. All this family of 
maps is called a construction of shift or simply a shift of hypersurface S along 
trajectories of dynamical system (1.1). 
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Definition 1.1. Shift f t : S — > St of hypersurface S along trajectories of Newtonian 
dynamical system with force field F is called a normal shift if hypersurfaces St 
are orthogonal to the trajectories of shift. 

In order to construct the normal shift we should, at least, choose initial velocities 
v(p) being perpendicular to initial hypersurface S, i. e. v(p) = v(p) ■ n(p), where 
n(p) is a unitary normal vector to S at the point p: 

(1.3) x fe | -x fc (p), z fe | = i/(p) -n fc (p). 

lt=o lt=o 

But this is not sufficient. Initial data (1.3) by themselves do not provide orthogo- 
nality of St and trajectories of shift for t ^ 0. We are to use other opportunities due 
to the choice of hypersurface S, choice of modulus of initial velocity v(p) on S, and 
choice of the force field F of dynamical system (1.1). If we choose hypersurface S, 
then, in the case of success, we would have the construction of normal shift realized 
only on some special hypersurface S (or in some special class of hypersurfaces). In 
paper [1] we left the choice of S to be arbitrary, and have concentrated efforts to 
the choice of function v{p) in (1.3). Then we found that the proper choice of v(p) 
on an arbitrary hypersurface S is possible only under some definite restrictions for 
the choice of force field of dynamical system (1.1). This became the origin for the 
theory of dynamical systems admitting the normal shift. It was developed in 
the series of papers [1-16]. On the base of these papers two theses were prepared: 
thesis for the degree of Doctor of Sciences in Russia [17] and thesis for the 
degree of Candidate of Sciences in Russia [18]. 

In [1~18] we restricted ourselves to the case of smooth hypersurfaces St with no 
singular points. However, in the process of shifting ft: S — > St some singular points 
can appear (they are called caustics). In particular, we can observe the collapse 
of St into a point at some instant of time t = to followed by a blow-up of this point 
into further series of smooth hypersurfaces for t > t$. Without loss of generality 
we can assume that to = 0. Then initial hypersurface S = {p } consisting of only 
one point p appears to be singular, and for t > we have blow-up of this point p 
into a series of smooth hypersurfaces S t . 

Definition 1.2. Blow-up f t : S — ► St of singular one-point hypersurface S = 
{po} along trajectories of Newtonian dynamical system with force field F is called 
a normal blow-up if for t > smooth hypersurfaces St are orthogonal to the 
trajectories of this blow-up. 

The idea to consider blow-ups of one-point sets in the framework of normal shift 
was suggested by A. V. Bolsinov and A. T. Fomenko when author was report- 
ing results of thesis [17] in a seminar at Moscow State University. Partially this 
idea was realized in [19]. In that paper was shown that Newtonian dynamical sys- 
tems admitting the normal shift of hypersurfaces are able to implement normal 
blow-up of any point p in Ricmannian manifold M. More completely the idea of 
A. V. Bolsinov and A. T. Fomenko can be realized in special investigation. This is 
the main goal of present paper. 

2. Geometry of normal blow-up. 

Let po be some point of Riemannian manifold. Let's consider normal blow-up 
of this point f t : {po} — > St along trajectories of Newtonian dynamical system 
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(1.1). Consider some particular trajectory of shift. For t = it passes through 

the point p . Let v(0) be the velocity vector 
corresponding to the time instant t = 0. If 
v(0) ^ 0, then this vector can be normalized 
to unit length: 



Unit vectors (2.1) for various trajectories be- 
long to the same tangent space T po (M). They 
can be identified with radius-vectors of points 
on unit sphere a in the space T po (M). 

Definition 2.1. The blow-up f t : {po} — > St 
of the po along trajectories of Newtonian dy- 
namical system (1.1) is called regular if velocity vectors v(0) at the point po are 
non-zero for all trajectories of this blow-up and if points corresponding to unit 
vectors (2.1) fill the whole surface of unit sphere a in T po (M). 

In the case of regular blow-up all hypersurfaces St possess spherical topology for 
sufficiently small values of parameter t ^ 0. Points q of the unit sphere a in T po (M) 
can be used to parameterize points of hypersurface St- In order to do it we shall 
write initial data determining regular blow-up of the point po as follows: 



(2.2) x k \ = x k (p ), x k \ = v{q)-n k {q) 

lt=o lt=o 

Trajectories of Newtonian dynamical system (1.1) fixed by initial data (2.2) de- 
termine a family of maps ft : <J — ► St being diffeomorphisms for sufficiently small 
values of parameter t ^ 0. 

Suppose that we have regular blow-up of the point po being normal in the sense of 
definition 1.2. Let's consider the solution of Cauchy problem (2.2) for the equations 
(1.1). This is the set of n functions x 1 (q,t), . . . , x n (q,t). Due to initial data (2.2) 
we can write Taylor expansions for these functions at the point t = 0: 



(2.3) x k {q, t) = x k (p Q ) + v{q) n k {q) ■ t + 0(t). 

Denote v(q) = u(q) ■ n(q). Vector v(q) E T po (M) has the meaning of initial ve- 
locity for trajectory that corresponds to the point q on unit sphere a. In terms of 
components of vector v(q) the expansions (2.3) can be rewritten as 

(2.4) x k (q,t)=x k (p Q )+v k (q)-t + 0(t). 

Let u 1 , ... , u n ~ x be local coordinates of the point q on unit sphere a in T po (M). 
Due to local diffeomorphisms of blow-up f t : a — > St they can be used as local coor- 
dinates on hypersurfaces St- Let's represent functions x k (q,t) and their expansions 
(2.4) in local coordinates u 1 , . . . , u™ _1 : 



(2.5) 



x k {u\ ... , u n -\ t) - x k ( Pa ) + v k {u\ ... , u 11 - 1 ) ■ t + 0{t). 
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Time derivatives of the functions (2.5) determine velocity vector on trajectories of 
blow-up, their derivatives in m 1 , ... , determine tangent vectors to S t : 

_ ^ dx k d _ ® x>t ® 

k=l k=l 

Change in any one of parameters u 1 , ... , u 11 ^ 1 leads to transfer from one trajec- 
tory of blow-up to another. Therefore vectors n, ... , r n _i are called vectors of 
variation of trajectories or simply vectors of variation. For the components 
of the vector v(f) and for the components of vectors Ti{t) from (2.5) we derive 

(2.7) v k {u\... ,u n -\t) = ^=v k {u\... ,u n - 1 ) + 0(l), 

(2.8) t 4 («,...,« ,t)_^-j._ i + 0(t). 

The normality condition for the blow-up ft : {po} St implies orthogonality of 
velocity vector v to all vectors T\, . . . , r n _i in (2.6). Let's write this condition as 
the condition of vanishing of scalar products (v | Tj): 

(2-9) (v|r^EE^(x 1 ,...,x")^g=0. 

k=l r=l 

Substituting the expansions (2.5), (2.7), and (2.8) into the equality (2.9), we deter- 
mine the asymptotics of left hand side of this equality as t — > 0: 

(2.10) £ $>fcr(po) v k (u\ . . . , u"" 1 ) ° { j • * + 0(t) = 0. 

fe=lr=l U 



Right hand side of (2.10) is identically zero. Therefore from (2.10) we get 
(2.11) £ E BkrtPo) v\u\..., u- 1 ) ^ j = o. 



du 1 

fe=l r=l 

Here gkr(po) = 9kr{x 1 {po) 1 ■ • ■ , £™(po)) are the components of metric tensor at the 
point po referred to local coordinates x 1 , ... , x n in M. Looking attentively at the 
left hand side of (2.11), we see that it is exactly the scalar product of the vector of 
initial velocity v(g) and the derivative of this vector with respect to parameter u l : 

„ 12) w „ K , W) .iS^C.^.a 

This equality (2.12) proves the following theorem. 

Theorem 2.1. for the regular blow-up ft : {po} - * of the point po along the tra- 
jectories of Newtonian dynamical system (1.1) to 6e normal it should be determined 
by initial data (2.2) with v(q) being constant: v{q) = const ^ 0. 
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Let's denote by vq the constant appeared in theorem 2.1. Then we can write the 
initial data (2.2) in the following form: 

(2.13) x k \ =x k ( P0 ), x k \ = v -n k {q) 

\t=o lt=o 

Further we state the following definition, which is central for the theory developed 
in present paper below. 

Definition 2.2. Newtonian dynamical system (1.1) with force field F on Riemann- 
ian manifold M is called a system admitting normal blow-up of points if for 

any point po G M and for arbitrary positive constant vq > initial data (2.13) 
determine the normal blow-up / t : p — ► S t of the point p a along trajectories of this 
dynamical system. 

Definition 2.2 was first formulated in paper [19]. Theorem 2.1 shows that the 
condition v(q) = = const built into definition 2.1 is absolutely inevitable. 

Let's substitute v{q) — v a into the expansions (2.5), (2.7), and (2.8) and take 
into account the above notation v(q) = v(q)-n(q). As a result we get the expansions 

(2.14) x k (u\... ,u n -\t)=x k (pv) + v n k {u\... ,u n ~ 1 ) ■ t + 0(t), 

(2.15) v k (u\... , u n -\t) = z/ n k (u\ . . . , u™" 1 ) + 0(1), 

(2.16) r k (u\ . . . , u-\ t) = ^^'^^ ■ t + 0(t), 

which hold due to initial data (2.13) determining blow-up of the point p along 
trajectories of dynamical system (1.1). 

3. Dynamical systems 
admitting normal blow-up of points. 

Definition 2.2 introduces new special class of Newtonian dynamical systems. 
According to the results of [19], it is not empty (see theorem 12.1 in [19]). In present 
paper we study this new class of dynamical systems introduced by definition 2.2. 

Let F be the force field of Newtonian dynamical system admitting normal blow- 
up of points. Then, according to definition 2.2, by choosing an arbitrary point 
Po G M and by fixing an arbitrary constant v$ > one can construct normal blow- 
up ft : {po} —* St- In order to study this blow-up we consider hypersurfaces S t with 
spherical topology, determine local coordinates u 1 , ... , u™ -1 transferred from unit 
sphere a to S t , and define vectors (2.6) on trajectories of this blow-up. Then we 
introduce the following scalar products: 

(3.1) p i = ( V |Ti). 

Such scalar products were already considered above in formula (2.9). In thesis [17] 
they were called the functions of deviation. Functions of deviations (3.1) are 
the measure of deviation of blow-up f t : {p } — ► St from normality. In the case of 
normal blow-up all these functions are identically zero: ipi = 0. 
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Vanishing of the functions of deviation ifi at the initial instant of time t = 
follows from initial conditions (2.13) regardless to the choice of force field F of 
Newtonian dynamical system (1.1): 

(3.2) ipA = lim^u 1 ,... ,u n ~ 1 ,t) =0. 

lt=o 

Indeed, as t — > vector of velocity tends to its limit value v(0) = v$ ■ n(q), while 
vector Tj tends to zero, this follows from the expansions (2.8) for its components. 
Hence scalar product ipi = (v | Tj) tends to zero. 

Apart from (3.2), identical vanishing of the functions of deviation in the case of 
normal blow-up implies vanishing of their time derivatives ipf. 

,„> . I difiju 1 ,... ,u n ~ 1 ,t) 

(3 - 3) H= = * lnn o oi = °" 

In calculating (pi for t ^ we can replace differentiation in t by the covariant 
differentiation with respect to parameter t along the trajectory of blow-up: 

(3.4) <pi = V m = V t (v | Tj) = (V t v | Tj) + (v | V t Tj). 

For Vtv we have V*v = F. This follows from the equations of Newtonian dynamics 
(1.1). Therefore formula (3.4) for (pi now is written as follows: 

(3.5) & = (F|Tj) + (v|V t Tj). 

Vector F has a finite limit as t — > 0, it is determined by initial conditions (2.13): 
F — > F(po,^o • n (<7))- While vector Tj tends to zero. Therefore first summand 
in (3.5) vanishes in the limit as t — > 0. Now consider vector VtTj in the second 
summand. Let's write components of this vector: 



f) T k n " 

(3.6) V t r? = ^ + EE r -^- 



Then let's calculate limit as t — > in formula (3.6), using the expansions (2.14), 
(2.15), and (2.16) for this purpose. As a result we get 



(3.7) h^^-Hsr — -■ 

Here n fe (u 1 , . . . , u n_1 ) are components of unitary vector n(q) being the radius vec- 
tor of the point q on unit sphere a in tangent space T po {M). Derivatives of the 
vector n(u 1 , . . . , u™ _1 ) in u 1 , ... , u™ -1 are coordinate tangent vectors to the sphere 
a in local coordinates u 1 , . . . , it" -1 . Let's denote these vectors by Ki, . . . , K„_i: 

fe=i 

Then the relationship (3.7) can be rewritten as 

(3.9) V t Tj -» vo ■ Ki(q) as t -f 0. 
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For vector of velocity v, as was mentioned above, we have the relationship 
(3.10) v -> u ■ n(g) as t -y 0. 

From (3.9) and (3.10) we derive vanishing of the second summand in right hand side 
of (3.5) as t — > 0. Indeed, vector Kj(g) is tangent to unit sphere a at the point q, 
while vector n(g) directed radially. These vectors are perpendicular to each other, 
their scalar product hence is zero. 

Thus, both summands in right hand side of formula (3.5) vanish as t — > 0, hence 
the relationship (3.3) holds. Similar to (3.2), this relationship is fulfilled due to 
initial data (2.13) regardless to the choice of force field F of the dynamical system 
(1.1). Therefore we consider analogous relationship for second order derivatives 

(3 - n) H= = £s w = °' 

which also should be fulfilled in the case of normal blow-up of point po- Let's 
differentiate the equality (3.5) with respect to t. This yields 

(3.12) ^ = (V t F | n) + 2 (F | V t Ti) + (v | V«Ti). 

Components of force vector F depend on double set of arguments: on coordinates 
x , . . . , x n of the point on trajectory and on components v , ... , v n of velocity 
vector of this point. This means that vector F depend on the point of tangent 
bundle TM. Such vectors are not embraced by the ordinary concept of vector field 
on a manifold. Therefore in paper [6] the concept of extended vector field was 
introduced. The concept of extended tensor field is its natural generalization. 

Definition 2.3. Extended tensor field X of the type (r, s) on the manifold M is 
a tensor- valued function that to each point q = (p, v) of tangent bundle TM puts 
into correspondence some tensor from tensor space TJ(p, M) at the point p on M. 

Smooth extended tensor fields constitute an algebra over the ring of smooth 
functions on tangent bundle. It was called the extended algebra of tensor fields 
on M. In the case of Riemannian manifold one can naturally define two covariant 
differentiations V and V in extended algebra of tensor fields on it. First was called 
spatial gradient, second was called velocity gradient. Covariant derivative 
V;F of the force vector in formula (3.12) can be expressed through corresponding 
gradients of extended vector field F. For the components of the vector VtF in 
formula (3.12) we have the following expression: 

n n 

(3.13) V t F k = J2^sF k v s + ^2v s F k F s . 

s=l s=l 

We shall not comment formula (3.13), and we shall not describe in details all things 
related with extended algebra of tensor fields (see Chapters II, III, and IV in thesis 
[17]). Technique of using extended tensor fields is assumed to be known to reader. 

Second covariant derivative V U Ti in formula (3.12) is expressed through and 
VtTj. The matter is that components of any vector of variation of trajectories 
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r in case of Newtonian dynamical systems satisfy the system of linear ordinary 
differential equations of the second order: 



n n n 



(3.14) m =l«=lr=l 



m—1 s—1 r—1 

n n 



S = l 8=1 

Taking into account (3.13) and (3.14), we can bring formula (3.12) to the form 



r=l V s=l / 

^ ( ^ s (V s F r + V r ig + fV s V s iU <. 



(3.15) . r=1 V 5=1 



r—1 \ s—1 s—1 / 

Here and everywhere below, aside with contravariant components of vectors, we 
use their covariant components obtained by lowering index by means of metric: 

n n 

Vi = ^2gijV 3 , F i = ^2 / g ij F : >. 

Quantities F^, VkF s , V s Fk in formula (3.15) are the components of smooth ex- 
tended tensor fields. They all have finite limits as t — * 0. Limits are determined 
by substituting local coordinates of the point po and components of the vector 
v(0) = Vo ■ n(q) for their arguments. Therefore in order to calculate limit of the 
derivative <pi it is sufficient to use the relationship (3.9) and remember that Tj — > 
as t — > (the latter is due to the expansions (2.16)): 



(3.16) Jim <fii = ^2vo 2 F r + ^ v s V r F s K\. 

r=l \ s=l / 

Substituting (3.16) into (3.11), we obtain the following relationship: 

n / n \ 

(3.17) Yl hF r + Y,v s V r F s \KT = 0. 

r=l V s=l / 

Here .KT are the components of the vector K; from (3.8). Note that left hand 
side of (3.17) is linear with respect to components of the vector Kj, while vectors 
Ki, . . . , K„_i form a base in the hyperplane perpendicular to the vector n(q). 
Vector n(q), in turn, is collinear to the velocity vector v(0) = v$ ■ n(q). Therefore 
if we introduce the operator P of orthogonal projection to the hyperplane perpen- 
dicular to velocity vector v and if we denote by PJ its components, we can replace 

(3.17) by an equivalent relationship 

n / n \ 

(3.18) J2\ 2F r + 12 vS VrF s )P[ = 0. 

r=l V s=l / 
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Orthogonal projectors P form an extended tensor field of the type (1, 1). Compo- 
nents of this field can be written in explicit form: 

(3.19) P[ = S r i -N r N i . 

Here <5[ is Kronecker delta symbol, while N r and Ni are contravariant and covariant 
components of extended vector field N formed by unitary vectors collinear to the 
vector of velocity v: 

(3.20) w=M, N = ^' 

What is the meaning of the derived relationships (3.18) ? The matter is that the 
relationships (3.11), in contrast to (3.2) and (3.3), cannot be fulfilled only due to 
initial conditions (2.13). They are equivalent to the relationships (3.18) that should 
be fulfilled at the point p for all vectors v such that |v| = v . If dynamical system 
(1.1) belongs to the class of systems admitting normal blow-up of points, as we 
assumed above in the beginning of this section, then the relationships (3.18) for its 
force field F are fulfilled at all points of tangent bundle TM, where |v| ^ 0. In this 
case they are partial differential equations with respect to the components of force 
vector F. 

Further we continue to study the relationships like (3.11). Next in the series of 
relationships (3.2), (3.3), and (3.11) is the vanishing condition for third derivatives 
of the functions of deviation ipi, ... , if n -i- 

(3.21) "ft! = fa =Q. 

\t=o «-o dt 3 

In order to calculate third derivative <p'j we differentiate the equality (3.15) with 
respect to t. The equality (3.15) has the following structure: 

n n 

(3.22) ^ = ^> r V t < + 5> r <. 

r—1 r—1 

Here a r and (3 r are components of extended covector fields. Therefore 

n n n 

Vi = J2 ar V « T ^ + 12 ( v * a - + VtT *' + 12 v *& T * = 

r—1 r—1 r—1 

n n / n n \ 

= a r Vttri + 12(12 Vs(x - yS + 12 ^ ps ) v *< + 

r—1 r—1 \ s—1 s—1 / 

n n / n n \ 

+ 12 & Vt< + 12 ( 12 V S + 12 H <. 

r—1 r—1 \ s—1 s—1 / 

In order to calculate Vta r and V^/3 r above we used formulas similar to (3.13). 
Further we take into account that r[(0) = (this follows from (2.16)). Then 

n n / n n \ 

"^(0) = a r v«<(0) + 12 [Pr + 12 v ^ yS + 12 ^ ps v *<(°)- 

r—1 r—1 \ s—1 s—1 / 
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Second covariant derivative V«r?" can be determined from the equation (3.14). In 
the limit as t — > this equation yields 

n 

(3.23) Vtt <(0) = £V s i^V t 7f(0). 

s=l 

Let's substitute (3.23) into the above expression for V'i(O). As a result we get 

n / n n n \ 

r— 1 \ s— 1 s— 1 s— 1 / 

The value of V^r 2 r for t — is determined from (3.9). Therefore the condition of 
vanishing of third derivatives (3.21) leads to the following relationship: 

n / n n n \ 

r—1 \ s—1 s—1 s—1 / 

Components of vectors Ki, ... , K„_i in the relationship just obtained can be 
replaced by components of orthogonal projector P. Arguments for doing this are 
the same as in replacing the relationship (3.17) by (3.18): 

n / n n n \ 

(3.24) ^ A + ^ V r F s a s + Y^ V s a r v s + ^ V s a r F s P[ = 0. 

r—1 \ s—1 s—1 s—1 / 

Now we are to substitute explicit expressions for a r and f3 r into the relationship 

(3.24) . They should be taken in comparing formulas (3.15) and (3.22) for (p\ 

n 

(3.25) 8=1 

But before doing this substitution, note that previously obtained equations (3.18) 
for the components of force vector F can be written as 

n 

(3.26) J2a r P[ = 0. 

r=l 

Let's apply the differentiations V and V to (3.26) and let's contract the resulting 
equalities with components of vectors v and F respectively This yields 



j2J2 vSv ^ p ^ = ' 



(3.27) r=1 8=1 

r—1 s—1 r—1 s—1 
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In deriving (3.27) we took into account (3.26) and we used the relationships 

(3.28) V.PT = 0, V S PT = ~ PI + J2 m PI N^j . 

The relationships (3.28) can be proved by direct calculations on the base of formulas 
(3.19) and (3.20) (see §5 in Chapter V of thesis [17]). 

Let's use the relationships (3.27) in order to simplify the equations (3.24). Due 
to the first of these relationships the third summand in (3.24) vanishes. Second 
relationship (3.27) enables us to transform fourth summand in (3.24). As a result 
of both these transformations we obtain 

(3.29) E & + E W a s + E )P[ = 0. 

r=l \ s=l s=l 1 1 / 

Now let's substitute a r and (3 r taken from (3.25) into the equations (3.29). Then 

n n n n 

E E V "^' F r + P l + E E F ^sFr PI + 

r— 1 s—1 r—1 s—1 

n n n n n 

(3.30) + E E 2 F * * r F' ^ + EEE w9 ^ s F q V r F° P[ + 

r—1 s—1 r—1 s—1 q—1 

71 n 9 TP ATS tp n n n 

r—1 s—1 ' ' r—1 s—1 q—1 

Further transformation of the obtained equations we begin with the fourth sum- 
mand in (3.30). Due to (3.19) we have <5[ = P[ + N r N { . Therefore 



n n n n n n n 



EEE^.^v r ^J7 = EEEE^.^iyv r jtfJ7^ 

r—1 s—1 q—1 r—1 s—1 q—1 j — 1 

n n n n n n n 

+EEEE^^^^^ FJ ^ = -EEE 2 ^^^ FJ ^'+ 

r—1 s—1 q—1 j — 1 r—1 s—1 j—1 

n n n n n n 

+Y,Y,Y,Y, Nq v s F q N s viv r F J p[ = -Y,Y, 2F ^rF s p:+ 

r—1 s—1 q—1 j—1 r—1 s—1 

n n n nun 

r—1 s—1 j—1 r—1 s—1 q—1 

Apart from the relationship S[ = P[ + N r Ni following from (3.19), here we used 
the relationship (3.26) written as 



E E M N ' VrFs P[ = ~J2 2F r PI- 

r—1 s—1 r—1 
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Let's apply this relationship once more for to transform second summand in the 
above expression. As a result we have 



n n n 



E E E vq ^sF q v r F s p[ = -J2Y, 2F * ^rF s pi 

r—l s—1 q—1 



r— 1 s—1 

n n n 



-EE lz ^^-EEE 2 ^ s w^. 

r— 1 s—1 r—l s—1 q—1 



In substituting this expression for the fourth summand into (3.30) we find that 
third summand cancels, while fifth and sixth summands change their signs. The 
equations (3.30) in whole now look like 



(3.31) 



r—l s—1 

n n 



r—l s—1 

n n n 



r—l s — 1 ' ' r—l s—1 q—1 



And finally, let's do some slight (purely cosmetic) transformations in the equations 
(3.18) and (3.31). Then write them combining into a system: 



(3.32) 



E :/ < • E v < (^F j ))pi = o, 



3=1 



E E + v ^ 2 v ~ 2 F * f j) Ni p * + 

i=l 3 = 1 



i=l 3 = 1 



r=l 



Let's state the result following from the above calculations in form of a theorem. 

Theorem 3.1. For Newtonian dynamical system (1.1) in Riemannian manifold M 
to admit normal blow-up of points its force field should satisfy the equations (3.32) 
at all points q = (p, v) of tangent bundle TM, where v ^ 0. 



4. Weak normality condition. 

Above we considered vanishing conditions for the functions of deviation ipi and 
for their derivatives ipi, (pi, and '<Pi at the initial instant of time t — 0. And we have 
found that first two conditions </?,(0) = and tfi(0) — are fulfilled only due to 
initial data (2.13) determining the blow-up of point. They make no restriction for 
the choice of force field F of dynamical system. Considering next two conditions 
<fi(0) = and '^(0) = 0, we derived the restrictions for F in form of the equations 
(3.32) for the components of force vector. 

Further we could step by step consider the vanishing conditions for the derivatives 
of the functions of deviation of higher order getting more and more equations for F 
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in each step. However, as we shall see soon, it is not necessary. The matter is that 
the equations (3.32) are exactly the same as weak normality equations, which 
arose in considering Newtonian dynamical systems admitting the normal shift of 
hypersurfaces. For the case M = M. 2 they were first derived in [1], then in [3] they 
were generalized for the case M = R n . In form of (3.32) corresponding to the 
case of arbitrary Riemannian manifold these equations were derived in [6] (see also 
Chapter V in thesis [17]). Weak normality equations are equivalent to the following 
condition of weak normality. 

Definition 4.1. Newtonian dynamical system (1.1) in Riemannian manifold M 
satisfies weak normality condition if for each its trajectory there exists some 
ordinary differential equation 



such that any function of deviation ip = (v | t) corresponding to the arbitrary choice 
of the vector of variation r on that trajectory is the solution of this equation. 

Words "any function of deviation" and "arbitrary choice of the vector of vari- 
ation" in this definition should be commented. Suppose that p = p{t) is some 
trajectory of Newtonian dynamical system (1.1). Let's include it into some (ar- 
bitrary) one-parametric family of trajectories p = p(u,t), so that for u = we 
would have p(0,t) = p{t). Parameter u can be introduced, for instance, by making 
dependent on u the initial data in Cauchy problem that fixes our trajectory p(t). 
In local coordinates the family of trajectories p = p(u, t) is given by the functions 



Derivatives of the functions (4.2) with respect to parameter u determine vector r(t) 
on the trajectory p(t): 



(compare with formula (2.6)). This vector is called the vector of variation 
of trajectories. This vector -r(t) (constructed as described above) is implied in 
definition 4.1. Its scalar product with the vector of velocity v is a function of 
deviation corresponding to it: <p = (v | r). 

It is easy to show that components of any vector of variation r{t) constructed as 
above satisfy the differential equations (3.14) (see paper [6] or Chapter V in thesis 
[17]). And conversely, any vector r(t) with components satisfying the equations 
(3.14) can be obtained by formula (4.3) in the above construction. Therefore words 
"arbitrary choice of the vector of variation" in definition 4.1 can be understood as 
the choice of an arbitrary solution of the system of linear homogeneous ordinary 
differential equations (3.14). Function of deviation 



(4.1) 



(p = A(t) ip + B(t) if 



x 1 = x x (u, t), 



(4.2) 



x n = x n (u,t). 



(4.3) 




n 



V (i)=(v|T)=x; 



Vk T 



.k 



fe=l 
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corresponding to this choice of r in the case of general position satisfies linear 
homogeneous ordinary differential equation of the order 2n (see theorem 6.1 in 
Chapter V of thesis [17]). In special cases (for special choice of force field F) the 
order of this equation can be lower. Definition 4.1 separates the case, when order 
of the equation for ip is 2. For this case in [6] the following proposition was proved 
(see also theorem 6.2 in Chapter V of thesis [17]). 

Theorem 4.1. Newtonian dynamical system in Riemannian manifold of the di- 
mension n 2 satisfies weak normality condition if and only if its force field satisfies 
the system of differential equations (3.32) at all points q = (p, v) of tangent bundle 
TM, where v ^ 0. 

Note that vectors of variation T\ 1 ... , t„_i arising in blow-up of points are 
naturally embedded into the above construction (4.3). Therefore if force field F of 
dynamical system (1.1) satisfies the equations (3.32), then corresponding functions 
of deviation ipi, ... , ip n -i satisfy the differential equations of the form (4.1). In 
this case the conditions 

(4.4) <pA =0, <pA =0 

lt=o lt=o 

from (3.2) and (3.3) provide identical vanishing of the functions of deviation. The 
conditions (4.4) by themselves, as we noted above, are provided only by initial data 
(2.13). Therefore we can strengthen the theorem 3.1 formulating it as follows. 

Theorem 4.2. Newtonian dynamical system (1.1) in Riemannian manifold M 
admits normal blow-up of points if and only if its force field satisfies the equations 
(3.32) at all points q = (p, v) of tangent bundle TM, where v^O. 

The restriction in dimension n ^ 2 from theorem 4.1 is inessential. We do not 
formulate it explicitly, since in the dimension n = 1 the concept of normal blow-up 
of points has no meaning. 

5. Concluding remarks. 

Theorem 4.2 is a central result of present paper. It reduces the study of Newto- 
nian dynamical systems admitting normal blow-up of points to the analysis of the 
system of partial differential equations (3.32) for their force fields. Moreover it pro- 
vides geometric interpretation for weak normality condition, reducing this rather 
abstract condition to visually obvious condition that dynamical system is able to 
perform normal blow-up of points. 

Note that in the dimension n = 2 (this case was studied in details in thesis 
[18]) the condition that Newtonian dynamical system admits the normal shift of 
hypersurfaces is also reduced to the system of the equations (3.32). Therefore we 
have the following proposition. 

Proposition 5.1. Class of Newtonian dynamical systems admitting the normal 
shift of hypersurfaces and class of Newtonian dynamical systems admitting normal 
blow-up of points for n — 2 do coincide. 

For the dimension n > 3 in the theory of Newtonian dynamical systems admitting 
the normal shift of hypersurfaces, apart from weak normality equations (3.32), we 
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have so called additional normality equations. They have the following form: 



i— 1 j—1 \ rn—1 / 

n n / n ~ 

i— 1 j—1 \ rn—1 

n n n n n r>j rpi ryrr 

EE^w = EEE "1^ ' 



Analysis of complete system of normality equations combined by (3.32) and (5.1) 
was undertook in [16]. However, in paper [16] an error was committed. Therefore 
part of results of [16] are not valid. This error was corrected in Chapter VII of 
thesis [17] (see also paper [19]). As a result an explicit formula for general solution 
of complete system of normality equations (3.32) and (5.1) was derived. 

Currently the analysis of separate system of weak normality equations (3.32) is 
urgent. In particular, would be worth to know whether something like proposi- 
tion 5.1 is valid in the dimension n ^ 3. Theorem 4.2 reduces this problem to the 
study of the equations (3.32) 
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